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Abstract 

We analyze the Poisson structure of the time-dependent mean-field 
equations for bosons and construct the Lie-Poisson bracket associated 
to these equations. The latter follow from the time-dependent vari- 
ational principle of Balian and Veneroni when a gaussian Ansatz is 
chosen for the density operator. We perform a stability analysis of 
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both the full and the linearized equations. We also search for the 
canonically conjugate variables. In certain cases, the evolution equa- 
tions can indeed be cast in a Hamiltonian form. 

PACS numbers : 05.30.-d, 05.30. Jp, 03.65.-w, 03.65.Fd 
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1 Introduction 



Poisson structures have already been used to study various properties 
of classical theories such as the Maxwell- Vlasov equations and the BBGKY 
hierarchy || [(J. They are useful in quantization methods; for a review, 
see @, |7|, |j. In the quantum case, besides the abstract quantum group 
approach, Poisson structures have found some applications in the many- 
body problem. In particular, they have been used to study the stability 
properties of important types of approximation such as the static Hartree- 
Fock (HF), the time- dependent Hartree-Fock (TDHF) and the random phase 
approximation (RPA), both in the zero temperature |J and in the thermal 
case [f|]. With their help, the connections between these theories, as well as 
with classical dynamics, follow in a very transparent fashion. 

The purpose of the present paper is to extend the analysis of ref.[|], 
which deals with fermions, to the time evolution of many-boson systems. 
Among the relevant problems that we have in mind, we can quote bose con- 
densation phenomenon [jTIJ, [TTJ, |T2[. In particular, recent experiments on the 
Bose-Einstein condensation of dilute atomic gases in laser traps at ultracold 
temperatures are raising a wide interest. They involve not only the 



study of equilibrium properties, but also that of dynamics. In the existing 
mean-field approaches, this dynamics deals with the expectation value of a 
creation or an annihilation operator, which is the order parameter, and more 
generally with the expectation values of product of two such operators. Our 
purpose is to analyze the structure of these equations so as to clarify their 
meaning. 
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Our paper is organized as follows. We briefly present in section 2.1 
the time-dependent variational principle proposed by Balian and Veneroni 
(BV)jl4[ [16] which we use to derive the approximate evolution equations. 
We also introduce, as well as some useful notations, the trial subspaces 
which we will consider throughout this work. These ansatze lead to the 
Time-Dependent Hartree-Bogoliubov (TDHB) equations. They belong to 
the class of the time-dependent mean-field approximations and they are the 
bosonic counterpart of the time-dependent Hartree-Fock (TDHF) equations 
for fermions. In section 3, we construct the Lie-Poisson bracket associated 
with these TDHB equations. This is achieved by extending the analysis of 
ref.Q to bose systems. Our formultation includes the possible existence of 
bose condensates. We derive in section 4 the corresponding linearized equa- 
tions and study their stability in the sense of Lyapunov. A question that 
arises naturally in this context is whether canonically conjugate variables 
can be found. We perform a canonical transformation allowing us to rewrite 
the TDHB equations in a form which, in some circumstances, is well suited 
to the discussion of this question. The last section contains some conclusions. 

2 The Bosonic TDHB Equations 



2.1. The time-dependent variational principle 

The BV variational principle fll , |i~6|| is devised to evaluate the expectation 
value Tr AD (t\) of a given observable A, measured at a final time t\, when 
the density operator Dq is known at the initial time to- To this aim, one 
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introduces the action-like functional 

'cU(t) 



I{A(t),V(t)} = TrA(t )D(t ) + f 1 dtTrV{t) 1^^- - i[A{t), H 

J to \ dt 



(2.1) 

where the variational quantities are the time-dependent operators T>(t) and 
A(t), the latter being subject to the final boundary condition 

A{t l )=A. (2.2) 

The symbol Tr stands for a trace over the complete Fock space (we take 
h=l). The desired quantity is then provided by the stationary value of the 



functional (|2.1|) , where H is the Hamiltonian of the system. 

When the allowed variations 5A(t) and 5V(t) are unrestricted, the sta- 
tionarity conditions of the BV functional (|2.1| ) lead to the exact Liouville- 
von- Neumann equation for T>(t) and to the (backward) Heisenberg equation 
for A(t). When restricted trial ansatze are chosen for A(t) and T>(t), the 
extremization of (|2.1|) yields approximate, and in general coupled, evolution 
equations for A(t) and T>(t). 
2.2. Trial spaces and the TDHB equations 

Let us first introduce the trial class of operators, which we shall use 
throughout this work for T>(t) and A(t). To this end, we define the 2n- 
component operator a 



/a 



a 



at 



(2.3) 



where a' and a denote boson creation and annihilation operators in a given 
single-particle space of dimension n. In terms of the a's, the usual commu- 
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tation relations can be expressed in the compact form (i, j = 1 . . . 2n) 

( l n \ 

[aj, acj] = T i:j , t = ia 2 = , (2.4) 

\ 1« J 

where l n is the (n x n) unit matrix and 02 the (2n x 2n) second Pauli matrix. 

With this notation, our trial class for the operator T>(t), the exponentials 
of the linear plus quadratic forms in a, can be written in the factorized form 



17, 19 



V = exp (z/) exp (Ara) exp (-arSa), (2.5) 

where v is a c-number, A a 2n-component vector and 5* a (2n x 2n) sym- 
plectic matrix (rS is symmetric). The tilda symbol denotes vector or matrix 
transposition. For any operator O, let us denote the average value (O) and 
the shifted operator 6 by (O) = TiOV/Z and 6 = O - (O). Then, V is 
completely specified by the knowledge of "the partition function" Z = TrD, 
the vector (a) and the contraction matrix p (which we shall also call the 
single-particle density matrix) defined as 

pij = ((Ta)jai). (2.6) 

For our purpose, the parameterization of the operator T> by 2, (a) and 
p is more convenient than that by u, A and S. In particular, we will see that 
the evolution equations for the former variables take a simple form. 

With the same notations, we take as our variational choice [T7], [Tj| for 
the operator A(t) a linear plus quadratic form in the operator a: 

A{t) = v a {t) + \ a {t)ra + -arS a (t)a, (2.7) 
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to which we shall loosely refer as a single-particle operator. 

With the choices Q2.5Q and (|2.7|) , it is possible to write explicitly the 
functional (|2.1|) in terms of the variational parameters Z, (a), p, v a) \ a and 



S a [fLY]]. Indeed, we find 

Ti AD(H) - £ it 

'§ 



J 



I . 1^4 - iZ \L a T{i{a) - t d£ 



m , -^S a (ip + 2\p, §]) 



(2.8) 



with L a = \ a — S a (a). The TDHB equations 

iZ = 0, 



i(a) 
ip 



T 



d£ 



d(a) 



(2.9) 



-2 



n d£ 



are then obtained [T7, IS] as the stationarity conditions of the reduced form 
( |2.8|) of the BV functional ( |2.1|) with respect to the parameters u a , X a and 
S a . In these equations, £ = (H) is the mean-field energy and the dots denote 
gj. We have used the usual convention (d/dp)ij = d/dpji. 

Thus, within the framework of the BV variational principle fl2.1|) , and 
within the trial class ( |2.5|) for the density operator, the time-dependent mean- 
field equations ( |2.9j ) appear to optimize the expectation values of single- 
particle observables. 

The mean-field energy is a sum of contributions of the type p and (a) (a) 
for its kinetic part, pp, p(a) (a) and (a) (a) (a) {a) for its potential part. Hence 



the equations (|2.9|) couple the dynamical variables (a) and p. If this coupling 
is disregarded, and if the interaction has zero range, the equation for (a) 
reduces to the Gross-Pitaevskii equation |T2]j. 



The equations ( |2.9| ) do not depend on the parameters Waif), A a (t) and 
S a (t): the evolution of the operator T>(t) decouples from the evolution of 
A{t). The knowledge of z/ a (t), \ a (t) and S a (t) is unnecessary for the deter- 
mination of T>(t), and therefore of TrAD(ti). The solution of eqs.( |2.9| ) is 
uniquely defined once the initial conditions Z(t ), (a)(t ) and p(t ) are spec- 
ified. It should be noticed that, in the present context, this feature results 
from the linearity in the variational parameters of the Ansatz ( [2.7] ). Other- 
wise, in the case of a non-linear parametrization, one should have to consider 
simultaneously the evolution equations for A(t) obtained by writing the sta- 
tionarity conditions of ( ^.1| ) with respect to the parameters characterizing 
T>(t). Then, the motions of T>(t) and A(t) would be coupled; moreover, the 
boundary conditions on T>(to) and A(t\) would lead to a mixed boundary 
value problem fT4], |18| . 

Although we shall not discuss them in the present paper, the equations of 
motion of A(t) have their own interest; indeed, they are tightly related to the 
random phase approximation (RPA) and they play a crucial role |2(| for the 
coherent determination of fluctuations and (two-time) correlation functions 
of the single-particle observables ( |2.7|) . 

Let us now discuss briefly some general propreties of the TDHB equations 



( |2.9|) (see [|T|, |17| for a more detailed discussion). We note first that the 
"partition function" Z is conserved; this is automatically entailed by the 
BV variational principle when the allowed variations 5A(t) includes the unit 
operator. Other properties of the TDHB equations are the conservation 
of the mean-field energy £ (for a time-independent Hamiltonian) and the 
unitary evolution of the contraction matrix p{t). The latter means that the 



eigenvalues of p(t) are conserved. This implies in particular the conservation 
of the single-particle von-Neumann entropy S = — TrPlogP and that an 
initially pure state p (p + 1) =0 remains pure during the TDHB evolution. 
The conservation of the free energy T = £—/3~ 1 S naturally follows is the 
temperature.) We will see the interest of this property in the next sections. 
These conservation laws show that the TDHB approximation reproduces, in 
the reduced single-particle space, some of the exact ones. 

3 The Lie-Poisson Structure of the TDHB 

Equations 

We now examine the mathematical structure of the TDHB equations 
^9f). A first interest of this analysis is to establish a close connection between 



the TDHB equations and dynamical systems [y, [21], |22[. It may prove useful 
in the elaboration of more sophisticated approximations. 

It was pointed out in ref.pi] that the TDHF equations possess a Lie- 
Poisson structure. The reason is that they deal with the expectation values 
of the one-body operators a\aj which form a Lie algebra. In the boson case, 
the set of the operators {1, a, aa} form also a Lie algebra. One expects, 
therefore, the existence of a Lie-Poisson bracket associated with the TDHB 



equations ( |2.9| ). This existence results also directly from a theorem which can 
be found in refs.jT], ^3|). In order to construct the corresponding Lie-Poisson 
bracket, let us introduce the set of the relevant operators 6 M (p = 1, 2, 3) 
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defined as follows: 



£>?„■ = (Ta)j(Xi . 



(3.1) 



6 1 = 1 , = 
These operators satisfy the Lie-algebra 

[6 m \?\=F?W. (3.2) 

Using the commutation relations ( j2.4| ), we find that the nonvanishing struc- 
ture constants f£ v are: 



<%J1 



iff 1 

^j'33\mn 



(fi 2 )ji,k = 8 ik (a)j - Tij(r(a)) k , 
{f m k 



(3.3) 



)i,jk 

12 >,.jk = (fDjiJ = SikSjl + TijTkh 

1 3 Jij,kl — ^il^jm^kn ^jk^lm^in ~l~ ^ik^jrnTln Tjl&inJkm- 

The Lie-Poisson bracket associated with the TDHB equations can now be 
constructed as follows. Consider the quantities z^ = which are directly 
related to our dynamical variables Z, (a) and p. The quantity z = (z l ,z 2 , z 3 ) 
is a vector in the Liouville space. This suggests defining the Poisson tensor 
{z^z"} as 

= {Z", z u } = f^z a . (3.4) 

Upon using ( |3.3| ), we can see that C^ u has only two non-zero components, 
namely 



■>22\ 



T iji 



(3.5) 



(C 33 )ij,ki= Supkj - 5 kj pu + Tik(rp)ij - Tij(pr)ik. 
Obviously, C^ v is antisymmetric with respect to its lower indices owing to 
the antisymmetry of the matrix r and to the property pji = (r(l + p)r)ij 
(seepBT). 
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We define the Lie-Poisson bracket of two functions h(z) and g(z), depend- 
ing on the components of the variables z^), in the usual way0: 

{M> = i^- (3-6) 

This expression is clearly bilinear and antisymmetric. In addition, it satisfies 
Leibniz's derivation rule 

{hg,k} = h{g,k} + {h,k}g, (3.7) 

and Jacobi's identity. This last property follows directly from the Lie algebra 
(|3.2|). Replacing in (RTBT) the function h by one component of the variables z 



and using (|3.5| ), we obtain for any function g 

{!,<?}= 0, 

\(^), (3.8) 

{pij,g}= -jip^h- 

In particular, when g is the energy or the free energy, we can express the 
TDHB equations in the form 

&* = = {z?,F(z)}, (3.9) 
which can also be written as 

z» = -C^^. (3.10) 

l oz 



The expressions ( |3.9| — pTT0| ) clearly exhibit the classical structure of the TDHB 



equations, with the mean- field energy £ , or the free energy J 7 , playing the role 
of an Hamiltonian. However, the non-symplectic Poisson tensor ( |3.5|) induces 
a non-canonical hamiltonian structure (also called a Lie-Poisson system). 
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From ( p.lO| ), one sees that the evolution of any different iable function 
h(z) is governed by the equation 

dh(z) 



dt 



{h{z),F(z)}. (3.11) 



In particular, any function h(z) which "Poisson-commutes" with the free 
energy (or the energy) is a constant of the TDHB motion. 

4 Linearized Equations and Stability 

Analysis 



Let z be a static solution of the TDHB equations ( |3.9|) : 



K,^(«>)} = 7 CMi >°)f| = o. (4.i) 

We consider a small amplitude motion 5z around z such that z^ = Zq + 5z^. 
Using ( |3.9| ) and retaining only linear terms in 8z, one obtains 



S^-T^W)" (4 ' 2) 

where JF 2 is the expansion up to second order of the free energy T: 

The linearized TDHB equations (LTDHB) can therefore be written in the 
form 

A (8z») = {5z»,F 2 (z)} , (4.4) 
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the Lie-Poisson bracket {, }o being given by ( |3.6| ) with the Poisson tensor 
evaluated at the point Zq. The LTDHB equations also exhibit a classical 
structure with JF 2 playing the role of an Hamiltonian (which we note is a 
constant of the LTDHB motion) and with the non-symplectic Poisson tensor 

The stability analysis of TDHB and LTDHB relies on the Lagrange- 
Dirichlet theorem (ref.[||] sect. 3.2). The TDHB motion is Lyapunov sta- 
ble (that is, a trajectory starting in the neighborhood of z$ remains in this 
neighborhood during the evolution) if Zq is a local minimum of T . The same 
condition (with T replaced by T<i) ensures the Lyapunov stability of the 
LTDHB motion provided the stability matrix Bn U = a u% v has only 

^ OZ OZ z=zq 

positive eigenvalues. In this case, the TDHB equations are said to be linearly 
stable. One notes that Lyapunov stability and linear stability are indepen- 
dent of each other since the non-linear terms in the TDHB equations may 
have stabilizing or destabilizing effects on the dynamics. 

Let us finally mention that the previous results hold both for pure states 
and at non-zero temperature. In the first case, one should evidently replace 
the Helmholtz free energy T by the mean-field energy £ . 



5 Search for Conjugate Variables 

We have seen in sections 3 and 4 that both TDHB and LTDHB admit a 
non-canonical hamiltonian structure. One may therefore wonder about the 
existence of canonically conjugate variables (CCV). This analysis is of course 
of great interest in the study of the classical limit or in the "re-quantization" 
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process of the mean-field approximation. 

In principle, the CCV are found by diagonalizing the Poisson tensor 
The second block C 22 being z- independent (see p75|) , one can easily find a 
global canonical parametrization for z 2 . Indeed, by using the usual cano- 
nical transformations (where j — 1 . . . n) 

aj = 71 ^ Xj + iPj " > ' fljt = V2 ^ Xj ~ iPi ^ ' ^'^ 

(which underlie the equivalence between the algebra of boson operators and 
the symplectic algebra of canonical pairs of the position and momentum 
operators x and p), and upon defining the variables 

* (5.2) 

TTi = (Pi), 

the evolution of (a) translates into the equations 

= dT 

~ ~d<p' 

which are clearly in an Hamiltonian form, and which, we recall, are coupled 
to p through T . 

In order to discuss the possibility of a canonical parametrization of z 3 , 
one must note that the block C 33 depends on the point z (more precisely, 
on p). One looks therefore for local CCV (there is no rigorous proof of non- 
existence, but the search for global sets of CCV seems very hazardous.) Let 
us consider the instructive (and simpler) example of the LTDHB equations 
in a basis where po is diagonal: poij = Poi^ij- The equation (|3.5| ) gives 
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and the LTDHB equations write 

i(Poi - Poj] 



dt 5pij 



dF 2 dT 2 N 

+ T - T 



d{5p) d(5p) 



ij 



It is easily seen that when p 0i — p j > 0, one can define the variables 

dpij + dpji _ i(S pij - 5 pji) 



Pij 



'2(poi-poj) ij2(poi-p j) 
of which the evolution takes the form of Hamilton's equations 

dF 2 . &F a 



(5.5) 



(5.6) 



Qii 



dpij ' '''' <)q;, ' 

poj (degeneracy in p), we see on 



However, for p Qi 



(5.7) 
that one must 



introduce new variables = Spij which are associated to the eigenvalue 
of the Poisson tensor. 

This scheme is general. When the Poisson tensor possesses vanishing 
eigenvalues (related to the existence of constants of the motion), it is un- 
avoidable, in order to obtain a canonical parametrization, to introduce addi- 
tional variables that Poisson-commute with the CCV. Let us illustrate this 
point by looking for a canonical parametrization of p in the pure state case 



(/T 1 = 0). Using 
matrices (see [|17| 



T]), the matrix p can be expressed in terms of the (n x n) 



A 



ij 



(PiPj), 

(Xipj+pjXi), 



(5.8) 



of which the evolution equations write 
f 



A 

A 



2(2fgr + A§£)„ 
-2 (2A|f + |f A), 



(5.9) 



- 2 [ a > i& 
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where the index s denotes the symmetric parts of the matrices involved. The 
pure state property p(p + 1) = is equivalent to the conditions 

hi =4(rA)i i -(A 2 ) ii = %, 
(Ar%, = (AI%, (5.10) 
(A A),, = (AA)y. 

As already mentioned in section 2, these conditions are preserved by the 
equations ( |5.9| ) if they were initially fulfilled. This is a typical example of 
vanishing eigenvalues of the Poisson tensor which depend on the initial con- 
ditions. Having chosen these initial conditions, one may in principle find a 
canonical form for the equations ( |5.9| ) with the additional constraints ( |5.10p . 
But in our particular case, we can explicitly eliminate these constraints by 
defining the two symmetric matrices 

f = r, 



G =|r- x A. 



(5.11) 



Indeed, the evolution of F and G is governed by Hamilton's equations 

p = d£ 

c " ^ 
G -—5F> 

and A will be given by A = |- F^ 1 + AGFG. It should be noticed however, 
that this canonical parametrization is not global since it does not span the 
whole manifold p, but is rather valid on a sheet characterized by the special 
initial conditions. 

Our last illustration in the search for CCV is the one-dimensional case. 
The matrices defined in (|5.8| ) become c-numbers. One can see that the quan- 
tity I is conserved, whatever the initial conditions. In particular, the initial 
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state needs not be pure. In order to rewrite the TDHB equations, it is possi- 
ble to keep the definitions (|5.11 ). Equations (|5.12p then remain valid (with 
A replaced by I/4F + 4FG 2 ). But sometimes, it is preferable to work with 
variables that have the dimensionalities of x and p. Hence, if we define 



w 
u 



we obtain from 



VT, 

1 A 



dT 
du ' 
_dT 
dw ' 



(5.13) 



(5.14) 



We arrive at the conclusion that a quantum one-dimensional bosonic system 
in the TDHB approximation evolves as a classical Hamiltonian system with 
two coupled degrees of freedom. In contrast with the n— dimensional case, 
this result holds for systems at non-zero temperature. 



6 Conclusions 



Upon selecting a trial time-dependent density operator T>(t) belonging to 
the class of gaussian bosonic operators and a trial time-dependent observable 
A(t) belonging to the class of "one-body" operators, we derived, from the BV 
variational principle, the evolution equations (labeled by TDHB) for the first 
three contractions of T>(t) namely, the "partition function", the expectation 
value (a) of the "one-boson" operator a and the contraction matrix p. 

We have investigated the relation of these TDHB equations with the 
hamitlonian formalism. We have exhibited the Lie-Poisson bracket associ- 
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ated with the TDHB equations and with their linearized form (LTDHB) by 
exploiting the Lie algebra of the set {1, a, aa} of the relevant observables. 
We have also performed a stability analysis of TDHB and LTDHB. 

The fact that the Poisson tensor associated with this bracket is non- 
symplectic, brings forward the main difference between the TDHB and Hamil- 
tonian evolutions, namely, that the former is non-canonical. This led us to 
wonder about the possible existence of canonically conjugate variables. The 
very existence of vanishing eigenvalues of the Poisson tensor, the dependence 
of the latter on the point of the manifold and, last but not least, the varia- 
tion of the multiplicity of the zero eigenvalues from point to point, render the 
search for global sets of CCV highly non-trivial in the general case. Local 
canonical parametrizations are however possible in spite of the additional 
constraints that one must introduce to account for the conservation equa- 
tions. This was shown to be the case for the LTDHB equations. We have 
also shown in this context, that the TDHB equations, in the pure state limit, 
can be cast in an Hamiltonian form by removing the constraints associated 
with the pure state conditions, and therefore reducing the number of degrees 
of freedom. 

The study of the mathematical properties of the TDHB approximation 
may also prove useful in the search for consistent extensions. For instance, 
one could "re-quantize" the TDHB equations by suitably deforming their 
associated Lie- Poisson bracket. 

We are indebted to R. Balian, M. Veneroni, C. Martin, H. Flocard and 
D. Vautherin for fruitful discussions. 
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